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Abstract 

We adapt continuous time random walk (CTRW) formalism to describe asset price 
evolution and discuss some of the problems that can be treated using this approach. 
We basically focus on two aspects: (i) the derivation of the price distribution from 
high-frequency data, and (ii) the inverse problem, obtaining information on the 
market microstructure as reflected by high-frequency data knowing only the daily 
volatility. We apply the formalism to financial data to show that the CTRW offers 
alternative tools to deal with several complex issues of financial markets. 
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1 Introduction 



The continuous time random walk (CTRW), first introduced by Montroll and 
Weiss (1965), has become a widely-used tool for studying the microstruc- 
ture of random process appearing in a large variety of physical phenomena. 
These range from transport in disordered media (e.g., Montroll and Shlesinger 
1984), earthquake modelling (e.g., Helmstetter and Sornette 2002) and even 
solar surface kinetics (e.g., Lawrence et al. 2001), to name just a few. This 
article extends the fields of application of the CTRW analysis by including the 
dynamics of financial markets. In this way, we aim to introduce the CTRW 
approach to the financial mathematics community where, to our knowledge, 
this formalism has been barely explored. 
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As its name suggests, the CTRW generalizes simple random walk models. Al- 
though the term "random walk" was coined by Pearson (1905), the formalism 
had been formulated in the seventeenth century in the context of gambling 
games such as the probability of ruin after betting n times in a coin tossing 
game. In this case, the sum of gains minus losses in n trials is equal to the 
state of a player's fortune (Weiss 1994). 

Financial markets have also been studied from the random walk point of view. 
In fact, this formalism was the first tentative model known in finance, hav- 
ing been suggested by Bachelier (1900) to describe stock market dynamics 
and give a price for a European call option. In fact, Bachelier modelled the 
price evolution assuming that prices change one unit at each time step with 
a probability p of going up and 1 — p of going down. Thus there are only two 
possible events. This process is called the binomial model and is the simplest 
random walk. Bachelier showed that the distribution of the resulting process, 
after a large number of time steps, tends to the Gaussian distribution. The 
limit behavior is often called the diffusion limit and the result obtained is just 
the consequence of the central limit theorem (Weiss). 

Several decades passed before there was further progress in the application 
of random walk methodology to analyze different aspects of financial transac- 
tions. Thus, in the book edited by Cootner (1964), there is a chapter devoted 
to the reexamination of the random walk hypothesis in finance. At that time, 
it was well-established that random walk models should be applied to the price 
return instead of the price itself as Bachelier asserted. In Cootner there are 
also several papers, basically published in the period 1961-1963, devoted to the 
question of whether the market is a random walk or a trend follower. Within 
these works we mention the article by Fama (1963) that studies an alterna- 
tive to the Gaussian Bachelier random walk, proposing instead the Pareto 
distribution. 

Later on, Cox and Ross (1976) used the Bachelier ordinary random walk 
formalism to provide a discrete time analog to the well-known Black-Scholes 
option price. In addition, Cox and Ross also obtained a different limiting 
process when the number n of time steps is large, finding that the binomial 
model leads to a Poisson jump process. Other contributions of the random 
walk approach to finance extend the binomial model by adding a third possible 
event, a crash, and observe the implications of it to the European option price 
(see e.g. Wilmott 1999). To our knowledge, there are not other financial market 
models exploiting all the possibilities that random walk analysis can offer in 
the study of many interesting phenomena in markets. 

Note that the ordinary random walk formalism mentioned above is based on 
the assumption that step changes are made at equal time intervals but this is 
a first approximation for many physical phenomena and markets. The CTRW 
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relaxes this restriction since it assumes that time intervals between transac- 
tions are not constant but random. Ticks have now, and in contrast with 
the Bachelier random walk, two sources of randomness: one coming from the 
amplitude and another one from the waiting times between ticks. This so- 
phistication is necessary to describe markets with tick-by-tick data (i.e., the 
highest frequency price data available). The deepest structure corresponds to 
the transaction-to-transaction operations. We believe that the CTRW formal- 
ism is especially useful in these cases and may improve the ordinary random 
walk approach. In contrast, it seems not to be suitable for cash indices since 
indices are averages of many prices and are recorded following some specific 
criteria described by an almost constant tick-by-tick frequency. 

Despite this promising fact, the CTRW is hardly known among financial an- 
alysts. Physicists have recently provided, from a new discipline called econo- 
physics, only a few examples of CTRWs applied to finance. Thus the papers 
by Scalas et al. (2000) and Mainardi et al. (2000) were among the first works 
addressed to this issue. Further developements were given by Raberto et al. 
(2002), Kutner and Switala (2003), and Masoliver et al. (2003). Perhaps one of 
the most solid reasons in favor of CTRW models is that they provide general 
expressions for the distribution of prices at time t in terms of two auxiliary 
densities that can be estimated from data: the probability density function 
(pdf) of the pausing time between ticks, ip(t), and the density for the magni- 
tude of the price increment at a given tick, h(x) (see Eqs. (J2]) and ([3]) below 
for a formal definitions of these quantities). 

Other quantities, such as the distribution of daily or longer-time prices based 
on two probability density functions ip(t) and h(x), can be obtained using 
the CTRW formalism. This in turn allows for the possibility of dealing with 
inverse problems, that is, estimating from the observed daily or longer-time 
data the forms of the microscopic functions ip(t) and h(x). This is useful since 
in many practical situations, one only has, at most, daily data. In this case 
our formalism enables us to determine features of the otherwise unknown 
microscopic structure of the financial process. 

Finally, we mention the existence of an important precedent that probably 
might have occurred to many readers. Within a very different approach, Mer- 
ton (1976) proposed a jump market model whose jumps are completely inde- 
pendent and have a Poisson distribution for the waiting times. The CTRW 
formalism generalizes the Merton model because we have now the freedom 
of choosing any probability distribution for pausing times and jumps and, 
even more importantly, the present formalism allows one to study arbitrary 
correlations between them. 

The paper is organized as follows. In Section [2] we present a short introduction 
to the CTRW formalism. In Section [3] we derive the exact distribution of 
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prices and volatility in terms of market microstructure statistics, while in 
Section H] we study the inverse problem. These sections contain analysis based 
on actual data to illustrate the procedures used. Finally, conclusions are drawn 
in Section and more technical questions are left to appendices available on 
the JEBO website. 



2 An introduction to the CTRW 

We first provide a definition of the CTRW from the perspective of financial 
markets. Let Z(t) be the return defined by Z(t) = \n[S (t) / S (t )} where S(t) 
is a speculative price and to is an initial time. In what follows we will assume 
that the return process is stationary so that Z(t) is independent of t . We are 
interested in the zero-mean return, X(t), rather than the return Z(t). This is 
defined by 

X(t) = Z(t)-(Z(t)), (1) 
where (Z(t)) is the average of Z(t). 

We now suppose that X(t) evolves following a CTRW and that any realiza- 
tion consists of a series of step functions (see Fig. [p. Therefore the return 
evolves discontinuously, and during any sojourn its value remains constant. 
In this picture X(t) changes at random times t , ti, t 2 , ■ ■ ■ , t n , ■ ■ ■. We assume 
that the intervals between successive steps, which we call sojourns or pausing 
times, T n = t n — t n _\ {n = 1, 2, 3, • • •) are independent, identically distributed, 
random variables with a probability density function, ip(t), defined by 

ip{t)dt = Prob{t < T n < t + dt}. (2) 

At the conclusion of a given sojourn X(t) undergoes a random change equal 
to AX n = X(t n ) — X{t n _i) whose probability density function is defined by 

h(x)dx = Prob{a; < AX n <x + dx}. (3) 

We collect these two random sources into one single function p(x, t) that rep- 
resents the joint probability density function of pausing times and return in- 
crements. That is, 

p(x, t)dxdt = Prob{a; < AX < x + dx; t < T < t + dt}. (4) 

In what follows, we take p(x, t) to be an even function of x thus assuring 
that there is no net drift in the evolution of X(t). Observe that two marginal 
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Fig. 1. Schematic representation of the return process. The dots mark the value 
X(t n ) of the return after each sojourn. T n = t n — t n -\ is the time increment of the 
n-th sojourn. 

densities can be formed out of p(x,t); the pausing-time density ip(t) and the 
density for a random jump h(x), these are related to p(x,t) by 



oo oo 

ipif) = J p(x,t)dx, h(x) = J p(x,t)dt. 

-oo 



(5) 



Our main goal is to obtain the probability density function of X(t). This 
function is the propagator, defined by 

p(x, t)dx = Prob{x < X(t) < x + dx}. 

It is well known that p(x, t), satisfies the following renewal equation (e.g., Cox 
1965 and Weiss 1994), 



t oo 



p(x, t) = po(x, t) + J dt' J p(x' , t')p(x — x',t — t')dx' 

-oo 



(6) 



where we have assumed that the initial jump occurred at t — 0. Equation 
involves several new definitions. The function p (x,t) is the propagator prior 
to the first jump and is equal to 



p (x,t) = *(t)6(x), 



(7) 



where is the probability that no transaction has occurred before time t, 



W(t) = J ^(t')dt'. 
t 
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The function S(x) appearing in Eq. (J7J) is the Dirac delta function that can be 
defined by 8(x) = for x ^ and 



5{x)f{x)dx = f(0) 



for any integrable function f(x) (Lighthill 1980). Equation ([6]) is derived from 
the consideration that at time t the first transaction has not yet occurred, this 
given by the first term on the right hand side of this equation, or else that a 
transaction occurred at time t' < t, at which time the return had value x' and 
from (x',t') the return process is renewed. We can solve Eq. in terms of 
the joint Fourier-Laplace transform: 

p(u,s)= I dte~ st I dxe luJx p{x,t) 



The solution is 

P{u,s) = - r, 9 

1 - p(u,s) 



where p (u;, s) and p(oo, s) are respectively the joint Fourier-Laplace transforms 
of the functions Po(x, t) and p(x, t). We easily see from Eq. (JTj) that the explicit 
form of po(u, s) is given by 

Po(^s)= 1 -^^, (10) 



where ip(s) is the Laplace transform of the pausing-time density ip{t). Equa- 
tion (J9j) furnishes the complete solution in the transform domain and is a 
convenient starting point in any CTRW analysis. 

One may argue that the CTRW is only applicable to intraday price evolution. 
The process can be extended to larger time lags as depicted in Fig. [2j Neglect- 
ing overnight effects on the price, obviously a first approximation, we show 
there how to describe the daily price evolution and thus represent the whole 
historical price evolution within the framework of the CTRW. This allows 
that the time t appearing in p(x, t) can be any time interval: minutes, hours, 
days, or longer. A significant effect of this linkage is that we unify within the 
same formalism high-frequency data with low-frequency data such as daily or 
weekly closing prices. We will see in future sections that this enables us to 
derive some important results. 
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Fig. 2. In (a) we describe the evolution of the stock price S(t). Consecutive transac- 
tions take place at discrete random times to, ti, ti, • • • and are represented by black 
circles. At every moment, the last traded price is the reference value during the in- 
tervals between changes. Two prices are noteworthy: the opening and closing prices. 
These are respectively the price of the first and last transaction and are represented 
by white (opening) and black (closing) boxes. In many situations opening and clos- 
ing prices are the only data available. Generally the closing price of a session does 
not coincide with the opening price of the next day. Nevertheless, as an approxima- 
tion, the last return of a given day can be taken as the first return of the next day. 
Therefore, since we do not model the behavior of the market when it is closed, we 
proceed as shown in (b); that is, we identify the last operation of a given day with 
the first transaction of the next day. In the figure both transactions are represented 
by white circles. 
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3 Prom market microstructure to return distributions 



We want to obtain the return distribution from information concerning tick- 
by-tick price behavior. Apparently, this is straightforward by taking the Fourier- 
Laplace inversion of Eq. (Q. Unfortunately, we need first to determine the form 
of the joint density p(x,t), but this is not possible from the available data. 
More easily accessible are the marginal pdf's, h(x) and ip(t), that fully define 
the tick-by-tick statistics. It is therefore essential to assume a functional re- 
lation between p(x,t) and its marginal densities ip(t) and h(x). Let us first 
study some possible ways of combining ip(t) and h(x) to obtain p(x,t) and 
discuss their main implications. 

The simplest choice would be based on the assumption that return increments 
and their duration time are independent random variables (Scalas et al.). In 
this case (ansatz 1) 

p(x,t) = h(x)i/>(t), (11) 



so that Eq. (jUJ) becomes 

,( U ,, )= N<*. (12, 
l-h(w)i/>(s) 



For the case of exponentially distributed pausing times, ip(t) = Aexp(— At), 
and after Laplace inversion, equation (IT21) results in the following characteristic 
function: 

p(u, t) = exp{-A[l - h(u)]t}. (13) 



However, in many situations one expects some degree of correlation between 
return increments and their duration, while Eq. (fTTj) implies complete indepen- 
dence between them. This is, for example, confirmed in the statistical analysis 
by Raberto et al. for the specific stock shares of General Electric. 

One possible assumption is to suppose that the density p(x, t) is such that its 
characteristic function p(u>,t) has the following functional form (ansatz 2), 



p{u,t) = if) 



h(u) 



(14) 



where is the pausing time density and h(u) is the Fourier transform 
of h(x). The assumption given by Eq. ( Tl4l) is at least intuitively plausible 
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because it implies that one must wait for a long time in order for a large 
variation of return to occur. In other words, large increments of the return are 
very infrequent. 

Let us show this by proving that sojourn times, T = t n — t n -i, and return 
quadratic increments, AX 2 = [X(t n ) — X(t n _i)] 2 , are positively correlated 
so that increasing return variations imply increasing sojourn times and vice 
versa. In effect, we define the following measure of the correlation between 
AX 2 and T by 



(AX 2 T) 
(AX 2 )(T) 



1. 



(15) 



We can easily evaluate the cross average (AX 2 T) using the joint characteristic 
function p(u;,t). Thus 



(AX 2 T) 



du 2 



tp(ui, t)dt 



LU = 



which after using Eq. (TU]) yields (AX 2 T) = 2(AX 2 )(T). Hence r = 1 as we 
meant to prove. 

We now come back to the ansatz in Eq. ( Tl4l) and observe that it allows us 
to write p(u, s) in the form p(u, s) = h(u)i/)[sh(u)]. This in turn leads us to 
write the formal solution to the problem given by Eq. ([9]) in the following more 
explicit form: 



p(uj,s) 



[l-${8)]/8 



1 — h(uj)ip sh{uj 



(16) 



We next consider an example in which it is possible to evaluate p(x, t) directly. 
Thus, on the assumption that ij)(t) = Aexp(— At) and h(x) = 7exp(— 7|a;|)/2, 
Eq. (USD yields 



X + s 



1 + 



A 



Xoo 2 /^ 2 + s 



(17) 



The inverse Fourier transform followed by the Laplace inversion and the con- 
volution theorem lead us to the following return pdf: 



p(x,t) 



S(x) + 4= f e?-^ 2 ^ 2 d£ 
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Fig. 3. The return pdf in dimensionless units p^x, Xt)/^ given in Eq. (| 18|) for 
x > 0, that is, without the delta function term appearing in (|18p . By the assumed 
symmetry in x, p(x, t) behaves similarly for negative returns. 

In Fig. [3] we plot this density for positive returns x > 0. 

The ansatz given in Eq. ffT^l) is not the only possible one that assumes a 
dependence between jumps and pausing times. Another plausible choice would 
be given by the following form (ansatz 3): 



Following the method outlined above one can easily see that in this case return 
variations and pausing times are also positively correlated but with a higher 
correlation (cf. Eq. ffT5]) ) given by r = 2. 

Note that in this case return pdf is given by 




(19) 




(20) 
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Fig. 4. The return pdf in dimensionless units p(jx, At) /j given in Eq. (|2ip for i/O, 
without the delta function term. 

and its joint Fourier-Laplace inversion yields 



p(x,t) = e- xt 5{x) + { 1 /2)e- xt/2 Q{Xt- 1 \x\)I U\H 2 - j 2 x 2 /2j , (21) 



where Q(x) is the Heaviside step function defined by 




when x > 0, 
when x < 0; 



and Io(x) is the modified Bessel function. We note that, contrary to the pdf's 
given in Eqs. ffl~3l) and ffl~8l) . the density p(x, t) given by Eq. fl2Tl) has a compact 
support since the value of the zero-mean return, X(t), differs from zero only 
within the interval —At/7 < x < At/7 for any fixed time t. Since this limitation 
on the value of X(t) is unrealistic, the model given by the ansatz in (jl~9l) 
together with the assumption that jumps are exponentially distributed seems 
not to describe events in real financial markets. In any case, for completeness, 
we plot this pdf in Fig. |U 
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3.1 Some asymptotic results 



We now focus on some asymptotic expressions of the distribution of returns, 
exemplified by its characteristic function p(cu,t). These approximate expres- 
sions are completely general and independent of the model chosen for p(x, t) 
as long as the sojourn time density ip(t) has some finite moments. Under 
these conditions, we have derived in Appendix A (see also Kotulski 1995 for 
an alternative method) some interesting results that can be summarized as 
follows: 

(a) If the jumps pdf h(x) has a finite second moment, fi2 = (AX 2 ) < oo, 
the asymptotic distribution of returns for long times approaches to the 
Gaussian density: 

p(u,t) ~ e'^ 2t/2{T) (t > (T». (22) 

(b) If h(x) is a long-tailed density, h(x) ~ as \x\ — > oo, then h(ui) 
has an expansion of the form h(u) — 1 — k\u\ a as uj — > (note that if 
< a < 2, h(x) has infinite variance). Moreover, if we assume that for 
uj small (Te lujAX ) ~ (T), then the asymptotic return pdf approaches the 
Levy distribution: 

p(u,t) ~ e - fc H^/m (23) 

(c) At intermediate times, t ~ (T), the behavior of p(x,t) for large values of 
|x| is the same as that of the jump distribution: 

p(x,t)~±h(x). (24) 

Note incidentally that the Levy distribution has been proposed by several 
authors as a good candidate for describing financial distributions when the 
observed fat tails are not Gaussian (e.g., Mandelbrot 1963, Fama 1963). How- 
ever, from the discussion above, we see that in the framework of CTRW's 
the appearance of Levy distributions is linked to the existence of transactions 
for which the jumps have infinite variance. Since this is manifestly impossi- 
ble in real markets, we conclude that the Levy distribution is an unsuitable 
candidate for the propagator. In order to overcome this inconsistency some 
authors propose a truncated Levy distribution (see, for instance, Mantegna 
and Stanley 1995), although in many cases the truncation procedure is done 
in an arbitrary way, without any justification provided by the behavior of the 
market. 
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3.2 The volatility 



Aside from the pdf p(x,t), which provides maximal information about the 
evolution of X(t), there is another quantity of considerable practical interest: 
the volatility, which does not require the knowledge of the entire jump dis- 
tribution h(x). It suffices to know the pdf ip{t) and the first two moments of 
h{x). 

Let us denote by rh 2 {s) the Laplace transform of the second moment (X 2 (t)): 



fh 2 (s) = / e~ st {X 2 (t))dt 



o 



This quantity can be derived in terms of the joint Fourier-Laplace transform 
of p(x, t) by 



m 2 [s) = - 



d p(u, s) 



duo 2 



(25) 



The combination of Eqs. (JH]) and (I2"5"j) leads to the relation 

*><•> = -i^frr, (») 

s[l-ip(s)\ 



where B,2(s) = — d 2 p(u , s) / dcu 2 



\lj=0- 



For the independent model defined by Eq. ffTTl) we see that B^is) = ^^(s), 
where fi 2 — (AX 2 ) is the second moment of the jump density. Hence 

m 2 {s) = ^—M— (27) 
s[l- ip(s)\ 



Note that the assumption of exponentially distributed sojourns, for which 
if)(s) = A/(A + s), leads to 

(X 2 (t)} = Xfi 2 t, (28) 



so that the variance always increases linearly with time as in the Wiener 
process. 

For the model specified by the ansatz we have R 2 (s) = fi 2 [sip(s)]' , where 
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the prime denotes a derivative, and 

M .) = ^f){ (29) 



For exponentially distributed pausing times we have m 2 (s) = /i 2 A 2 /s 2 (A + s), 
so that 



Xt + e- xt - 1 



(30) 



The limiting behaviors of this volatility are (X 2 (t)) ~ /i 2 (At) 2 for At <C 1, and 
(X 2 (t)) ~ ii 2 Xt for At > 1. 

For the second ansatz given by Eq. (}T9]) we see that ^(s) = /^^(s) 3 an d 
s 1 - 



In the case of exponentially distributed pausing times we have 

m 2 (s) = A 3 /i 2 /s 2 (A + s) 2 

and 



(X 2 (t))=to [Xt (l + e ~ xt )-2(l 



-xt 



(32) 



In this case the limiting behaviors are given by (X 2 (t)) ~ /z 2 (Ai) 3 /6 for Xt ^1 
and (X 2 (t)) ~ /x 2 At for At ^> 1. Therefore, and similar to ansatz ( Tl4l) . if we 
suppose that sojourn times are exponentially distributed, the volatility has 
a superdiffusion-like behavior at short times and a diffusion-like behavior at 
long times. 

We remark that the asymptotic variance grows linearly with time regardless 
of the specific form of p(x,t). Indeed, assuming that in all realistic models fi 2 
and (T) are finite from Eq. (1221) we readily obtain 

{X\t))^^-t (t»(T)), (33) 



independent of the correlation model. 
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3.3 The U.S. dollar/Deutsche mark case 



Let us now show how to infer the return statistics at an arbitrary time t from 
high frequency data. The case we discuss is that of the U.S. dollar/Deutsche 
mark future exchange, based on a transactions data base from January 1993 to 
December 1997. This consists of 1,048,590 transactions. We have chosen this 
data since this future has a high trading volume that is required to provide 
reliable statistical inferences. 

The first assignment is to estimate from high-frequency data plausible forms 
for ip(t) and h(x). In Fig. Owe plot the experimental pausing-time density 
tp(t) for the U.S. dollar /Deutsche mark future. We see that the data is well fit 
by a power-law density: 



where a = 3.47 and A = 2.73 x 10 2 s 1 . The mean sojourn time and its second 
moment are 

(T) = ^ (T 2 ) = (35) 

a — I (a — 2){a — 6) 



For the dollar/mark future market the experimental mean sojourn time di- 
rectly evaluated from data is (T) exp = 23.6 s, in satisfactory agreement with 
the theoretical prediction of (T) = 24.9 s evaluated from Eq. fl35|) with the 
parameters estimated from the fit to ip(t) in Fig. [51 

In Fig. [6] we plot the experimental jump density, h(x). The data indicate that 
the experimental h(x) can be considered a symmetric function of the return 
increment x. A good fit is also given by a power law with a greater exponent 
than that of ijj(t): 

KX) = 2 7 (1+~M/7F' (36) 



where (3 = 5.52 and 7 = 2.64 x 10~ 4 . Power law densities such as ( 1341) and (1361) 
have been suggested for describing several market models such as individual 
companies by Plerou et al. (1999), market indices by Gopikrishnan et al. (1999) 
and Mainardi et al. (2000). 

We have now all the elements to get the return pdf p(x, t) at any time t and, 
hence, to obtain a complete analysis of the dollar/mark exchange market. In 
order to focus ideas and results, we will now assume that pausing times and 
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Fig. 5. Empirical distribution of the time between the closest-to-maturity transac- 
tions of the U.S. dollar/Deutsche mark future (in the US market). The analyzed 
data range from January, 1993, to December, 1997. The pdf of the sojourn times, 
if)(t), clearly follows a power law at sufficiently long times. The solid line is the 
representation of Eq. (|34p . The normalized time between sojourns is t/(T). 

jump sizes are positively correlated as specified in Eq. (TT4"|) . This is supported 
by empirical data since for the U.S. do liar /Deutsche mark future market we 
have a positive correlation given by r = 0.53 (cf. Eq. (|T5|) ). 

Unfortunately, the densities given by Eqs. (|3^|) and (1361 do not allow an exact 
evaluation of p(x, t), or even of the characteristic function p(u, t). The Fourier 
and Laplace transforms of h(x) and ip(t) are too complex to manage them 
analytically as required for deriving the characteristic function ( JT6l) . We will 
therefore restrict ourselves to study the behavior of the tails of the distribution. 
We have shown above that, for moderate times, the tails of p(x, t) are the same 
as that of the jump density h(x). That is, for large |x| (cf. Eq. (EE 



p(x,t) 



(T) 



h(x) 



(*«<r» 



In our case we see from Eqs. (15BT)- (|5oT) that (T) = A 1 /(a — 2) and h(x) 
{(3 — l)7 /3_1 |x| _/3 /2 as \x\ — * oo. Therefore, 



p{x,t) 



(g ~ l)At t_ 
2(a-2) \x\ 



(3-1 



(t « <T» . 



(37) 



Hence, the tails of p(x, t) follow the same power law as that of the return 
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Fig. 6. Empirical distribution of the logarithmic changes between transactions in 
the U.S. dollar/Deutsche mark future market. Positive variations (increments) and 
negative variations (decrements) exhibit approximately the same behavior, thereby 
supporting our assumption of the symmetry of h{x). The plot suggests a power law 
dependence at large return increments. It also includes a graph showing the shape 
of h(x) in Eq. (|36j) . using the parameters reported there. 



increment distribution h(x). Observe that the approximation ( l37j) is valid 
when the time lag is of the same order of the mean sojourn. This prediction 
of the theoretical model is confirmed by actual data. In Fig. [7] we show the 
empirical p(x,t) for the dollar/mark future exchange for t — 15 seconds in 
which the time lag t is of the same order than the empirical mean sojourn 
(T) e xp = 23.9 s. The empirical distribution clearly shows a power-law decay 
with exponent (3 ~ 5.5, which coincides with the decaying exponent of h{x) 
which is consistent with the predictions of the CTRW model. 

We close this section with a remark about the volatility. In this case, we need to 
obtain the Laplace transform of the expression for ip(t) given by Eq. fl34|) . This 
results in an expression in terms of the Kummer function U(a, c, z) (Magnus et 
al. 1966), which renders very difficult the inversion of ^(s) to obtain (X 2 (t)). 
We thus limit ourselves to another asymptotic result. As we have proven above, 
the variance grows linearly with t as t —>■ oo regardless of the model chosen; 
therefore, {X 2 (t)} ~ /i 2 t/(T) (see Eq. (155])). In our case (T) = X~ 1 /(a - 2) 
and (X 2 (t)} ~ (a - 2)^ 2 At for At > 1. We refer the reader to Fig. HOI of Sect. 
14.11 where this asymptotic linear dependence on time is clearly observed for 
the U.S. dollar /Deutsche mark futures market. 
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Q. 




X 



Fig. 7. The empirical probability density function p(x, t) for a time lag t of 15 seconds 
(U.S. dollar /Deutsche mark future market). The model leads to a power-law decay, 
governed by an exponent (3 that is precisely the one appearing in the power-law for 
the jump density h(x) as expected when t « (T) exp = 23.9 s. We show the tail only 
for positive increments; the tail for negative increments behaves similarly. 

4 The inverse problem 

To this point we have determined the return density p(x, t) when the densities 
ip(t) and h(x) are known. However, from the CTRW formalism, it is possible 
to infer the mechanism of price formation at a microscopic level. This requires 
one to obtain the statistics of high frequency data represented by the pausing- 
time density ip(t) and the return increment density h(x) once we know the 
variance (X 2 (t)) and the return pdf p(x, t). We call this "the inverse problem" . 
Obviously this problem will require the knowledge of the complete pdf p(x, t) 
for all times t, something that is beyond reach, since, in practice, one has 
access, for instance, to daily data from which one is able to guess the pdf for 
the daily return p(x,t), where t = 1 day. In such a case we do not have an 
analytical expression for p(x, t) permitting us to get the microscopic densities 
ip(t) and h(x) exactly. Nevertheless with the formalism presented below we 
will be able to guess which forms of ip(t) and h(x) may be consistent with the 
observed p(x, t). 

Let us first discuss the problem of obtaining the microscopic jump density 
once we know the return pdf p(x,t). In effect, setting s = in Eq. (EJ) and 
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Fig. 8. The inverse Laplace transform of the pausing time distribution (|4ip . (|49p 
and (|54f) with r defined by Eq. (|42f) . The values of alpha chosen are the ones obtained 
from IBM (a = 0.915), 3M (a = 0.877), Gillette (a = 1.12) and Dupont (a = 1.06) 
data (see Tabled]). 

taking into account that p(u>, 0) = h(u>) and Pq(uj, 0) = (T), we have 



p(w,0) 



(T> 



1 - /i(u;) 



so that 

= 1 - 



<T> 



p(",0)' 



(38) 



Since 

oo 

p(u,0) = J p(u,t)dt, 



we see that Eq. fl38|) allows us to get the characteristic function corresponding 
to the jump density h(x) if we know the entire return pdf p(x, t) for all t > 0. 
Therefore, Eq. (1551) is very difficult to implement because, as mentioned, we 
know at most a numerical series for the value of the propagator instead of an 
analytical expression of p(x, t). In these cases, Eq. (|38|) is useless for evaluating 
h(x). However, Eq. (138]) can be useful for testing any statistical hypothesis and 
evaluating statistical parameters for the form of h(x). Incidentally, we note 
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that Eq. ([58]) is valid independent of the existing correlations between jumps 
and pausing times, that is, regardless of the form of the joint density p(x,t). 



Unfortunately, the inverse problem for obtaining the pausing time density ip(t) 
from the return variance (X 2 (t)) requires the knowledge of the form of p(x, t). 
Thus, for the independent model in which p(x,t) = h(x)ijj(t) the volatility is 
given by Eq. ( 1271) . for which we get 

- , s srhois) 



As an illustration of this procedure let us see for which density ip(t) corre- 
sponds to an anomalous diffusion volatility of the form 

(X 2 (t))=kt a (k,a>0) (40) 
In this case, 

m 2 (s) = kT(l + 
where T(x) is the Gamma function. From Eq. fl39l) . we get 

^ = TT aim I V (41) 

1 + p 2 s a /kT(l + a) 

We can invert numerically this Laplace transform. To this end we use an algo- 
rithm proposed by Stehfest (1970) that is the simplest method for numerically 
inverting Laplace transforms in the literature. We show in Fig. [8] the paus- 
ing time distribution for several values of a in terms of a normalized time r 
defined by 



T 



[kT(l + a)/p 2 ] 1/a t. (42) 



One can also obtain, by using Tauberian theorems (Handelsman and Lew 
1974), the asymptotic behavior of ip(t) as t — > oo and also the short time 
behavior t — > 0. It is worth to mention that Kotulski gives an alternative 
method based on renewal theory and limit theorems for random sum of jumps. 
In any case, following Tauberian theorems, we see from Eq. (l4T|) that for small 
s, 

kl [1 + a) 

and the Tauberian theorems imply that 

Ht) ~ / ' 2Si , n(7m) r 1 -" (kt a » M2 /r(l + a)) . (43) 

K7C 
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On the other hand, for long values of s, we have 

4>(s) ~ fer(1 + a) s~ a (s a > £T(1 + a)//i 2 ) . 

^2 



Hence, 



^( t )~^ t -i+« (H Q </i 2 /r(l + a)). (44) 



We thus see that a power-law volatility leads to a power-law pausing density 
at short and long times. 

Let us treat now the cases in which there exists a correlation between paus- 
ing times and return increments. When the ansatz in Eq. ffMl) is valid, the 
formulation of the problem given by Eq. ([TBI) provides a relation between the 
Fourier and Laplace transforms of the daily return density p(x, t) and the high 
frequency densities ip(t) and h(x). A second relation is required to determine 
the microscopic densities. This is supplied by the volatility. In effect, we start 
from Eq. (|29|) . which, on defining 0(s) = sip(s), can be written in the form 

d<p(s) „ . . 7, . „ , . 

ill — : hm 2 sW) S = sm 2 (s). 

as 

The solution to this equation with the initial condition 0(0) = yields 

${s) = —[ ^2(0 exp -- / m 2 (£)d£ de- (45) 



Equation ([431) yields an expression for the Laplace transform of the transac- 
tions pausing-time density once we know the volatility. 

In the case of a linear variance of the form 

(X 2 (t))=aH, (46) 

for which 77^2(5) = cr 2 /s 2 , we see from Eq. fj4"5]) that 

$( s ) = -—ePt'Eii-D/sY 
s 

where D = a 2 / n 2 , and Ei(z) is the exponential integral. The Laplace inversion 
is (Roberts and Kaufman 1966) 



i){t) = 2DK (2VDt), (47) 
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where Kq(x) is a modified Bessel function of the second kind (Erdelyi 1954). 
The asymptotic form of if)(t) is (Magnus et al.) 



i>{t) ~ — — e- 2vm (Dt > 1). (4* 



Therefore, a linear variance implies that ln[^(t)] falls off assyptotically as 
-VADi. 

Another example is provided by the anomalous diffusion volatility given in 
Eq. (jlD]). Now from Eq. (H5J) we get 

\kT(a)] 1/a 

^( s ) = i/a U[l/a, 1/a, kT(a)/^ 2 s% (49) 

M 8 



where U(a, c,x) is a Kummer function (Magnus et al.). This function can be 
inverted numerically with results shown in Fig. [SJ As before we use Taube- 
rian theorems to obtain the asymptotic behavior of ip(t). We have shown in 
Appendix B that 

Ah 

^{t)~ — t a ~ l (kt a <n 2 ), (50) 

fJ>2 



and 



KIT 



(51) 



where A is a constant given by Eq. ( IB. 21) . Note that these results agree with 
Eqs. (1431) and 

Let us briefly comment on the inverse problem under the assumption that 
pausing times and return jumps are correlated as in ansatz (fl9l) . In this case 
from Eq. (|26|) we have 



771,2 (S) 



\-${s) 



(52) 



Consequently, if we know the return moment (X 2 (t)) and hence its Laplace 
transform 7712 (s), the Laplace transform of the microscopic pausing-time den- 
sity ip{t) obeys the cubic algebraic equation 

l^2w{s) + sm 2 (s)ip(s) - srh 2 (s) = 0, (53) 
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whose solution is 



fa) =P(s) 1/3 \ 1 + ^1 + 8/3(^/27 



1/3 



+ 



1- Vl + 8/3(s)/27 



1/3' 



(54) 



where /3(s) = sm 2 (s)/2/i 2 . The numerical inversion of this equation when 
(X 2 (t)) = fct Q is shown in Fig. EJ A gain we can proceed as before and perform 
an asymptotic analysis based on the use of Tauberian theorems. We present 
the results of these calculations without going into details. We have 



kT(l + a) 
/x 2 r(a/3) 



r 1+Q / 3 {kt a < /i 2 ) 



(55) 



and 



/i 2 sm(7raj X _ Q 



(fct a >/i 2 ), 



(56) 



which is similar to Eqs. (J4"3l) - (j4"4"j) and fl50l) - fl5Tl) . We therefore see that the ap- 
pearance of asymptotic power-law pausing densities seems to be independent 
of the existence of correlations between jumps and pausing times. We also ob- 
serve the asymptotic distribution at long times has the same decay exponent 
for dependent and independent models. 



4-1 A market application 



We finally present an example involving real data. We base our discussion on 
four sets of data points corresponding to daily opening and closing prices for 
shares of IBM, Gillette, Dupont and 3M. These were chosen because a con- 
siderable amount of data was available for them. Thus data for IBM consists 
of daily prices from Jan. 2, 1968 to Dec. 31, 2002 with a total of 8,784 data 
points. The data for Gillette corresponds to daily prices from July 27, 1980 to 
Dec. 31, 2002 with 5,628 data points. Analogously the data for Dupont and 
3M consist of daily prices from July 28, 1980 to Oct. 17, 2003 with respectively 
a total of 5,801 and 5,842 points. For all these firms we have tried to deduce 
the form of the microscopic pausing time density, ip(t), horn a knowledge of 
low-frequency volatilities. We first plot (X 2 (t)) for t > 1 day. This is shown 
in Fig. [9j In the case of all of these companies we see a slightly but clear 
anomalous diffusive behavior for few days and a regular diffusive behavior at 
long times: 

(X 2 (t))~t a (t^0) and (X 2 (t))~t (t->oo), (57) 
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Fig. 9. The return variance dependence with time for IBM, Gillette, 3M, and 
Dupont. 



where a = 0.915 ± 0.005 (IBM), a = 1.12 ± 0.04 (Gillette), a = 1.06 ± 0.03 
(Dupont) and a = 0.877 ± 0.006 (3M). We therefore see that the anomalous 
diffusive behavior in the data is subdiffusive for IBM and 3M and superdiffu- 
sive for Gillette and Dupont. We summarize all the estimations of these fits 
in Table [U 

What can we say about the microscopic density ^(t)? From the analysis per- 
formed above we see that, in all cases, ip(t) dacays exponentially at long times. 
For intermediate times, times not too long neither too short, for instance 
t > {^/k) 1 ' 01 (see Eq. fj4"0]) ). we cannot apply either Tauberian or Abelian 
theorems. In these cases, as we have done above (see Fig. [8]), we must numer- 
ically invert ip(s) to obtain if)(t). 

Let us remark that a good empirical test of the formalism used for the inverse 
problem is to estimate the waiting time density. Unfortunately, we could not 
get enough data to perform such a test for the four stocks studied above. 
However, we can take the U.S. dollar/Deutsche mark high frequency data and 
filter them in order to get daily return changes. We thus observe the time- 
dependence of the variance provided by this daily data and perform a fit to 
the form kt a . Figure fTUI plots the resulting variance and its fit. 
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Table 1 

Summary of the linear and polynomial regressions for data from IBM, Gillette, 
3M and Dupont. The details of the linear and polynomial regressions are shown 
specifying the range within each fit is performed and the values with their errors 
obtained from each regression. Some parameters are given in days (d). 





IBM 


Gillette 


3M 


Dupont 


Period 


1968-2002 


1980-2002 


1980-2003 


1980-2003 


Points 


8,784 


D,OZO 


0,o4Z 


D,oUl 


Linear Regression: (X 2 (t)) = At + B 








Range in days 


15-25 


1 1 90 


1 ^ 9^ 
10— ZD 


1 ^ 9^ 
10— ZO 


(A ± 5 A) x 10 4 d" 


1 2.5 ±0.7 


8.104 ±0.005 


1.42 ±0.03 


2.17 ±0.01 


(B ± 5B) x 10 4 


7.30 ±0.14 


4.15 ±0.07 


7.1 ±0.5 


4.6 ±0.2 


Polynomial Regression: (X 2 (t)} = kt a 








Range in days 


1-4 


1-4 


1-5 


1-5 


(k ± 5k) x 10 4 d" 


- 1/a 2.129 ±0.010 


5.8 ±0.2 


2.73 ± 0.02 


2.61 ±0.07 


a ± 5a 


0.915 ±0.005 


1.12 ±0.04 


0.877 ± 0.006 


1.06 ±0.03 



Once we have estimated the a, we can numerically invert ^(s) and finally we 
can compare the pausing time distribution provided by the inverse method 
with the one directly observed from high-frequency data. We compare them 
in Fig. [TT] showing a good agreement in the region of intermediate times for 
which kt a > fi2- 

If we take a closer look at the three pausing distributions provided by the 
independent model and by the dependent models exemplified by ansatzes 1 
and 2, we observe that both the independent model and ansatz 1 describe 
the U.S. dollar/Deutsche mark data much better than does ansatz 2. This is 
explained by considering the quantity r defined in Eq. (1151) that quantifies 
the correlation between waiting times and jumps. We have seen that for the 
independent model r = 0, ansatz 1 gives r = 1, and ansatz 2, r = 2. On the 
other hand, the empirical correlation of the US dollar /Deutsche mark data 
is (see Section I3.3j) r = 0.53, and intermediate value between those of the 
independent model and ansatz 1, while the correlation given by ansatz 2 is 
higher. All of this seems to indicate that the actual Foreign Exchange (FX) 
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Fig. 10. The U.S. dollar /Deutsche mark daily return variance ranging from 1 day 
to 25 days return changes. We have filtered the tick-by-tick futures data to provide 
the daily return changes from January, 1993, to December, 1997, 1,035 days. Again 
the short time scale is fitted assuming that the variance has the time dependence 
kt a . The polynomial regression is done with nine data points (from 1 day until 9 
days) and yields k = (2.09 ± 0.05) x lO^days) 1 /" and a = 0.94 ± 0.03. The long 
time fit assumes a linear behavior, that is, At + B. The linear regression takes eleven 
data points (from 15 to 25 days) and yields A = (1.42 ± 0.03) x 10 _4 (days) _1 and 
B = (7.1 ±0.5) x 10" 4 . 

market should be described with an intermediate model between the indepen- 
dent model and ansatz 1 while the model based on ansatz 2 is too correlated. 
This is exactly what is observed in Fig. [TTJ since the curves corresponding to 
the independent model and ansatz 1 oscillate very closely around the empirical 
distribution from r = 0.1 to r = 20 while the curve corresponding to ansatz 
2 clearly deviates from the empirical result. 

We can quantify this by performing the chi-square test on these fits. The 
results, based on Fig. [TTJ, are the following. The value of chi-square statistic for 
the independent model is 6.5 times higher than that of the fitted distribution 
given by Eq. ([Ml) . For ansatz 1 the ratio is 39.2, and for the ansatz 2 the value 
is 60.7 times higher than that of the pausing fit. This is quite consistent with 
the assertions in the previous paragraph. 

We finally observe from the inset of Fig. [TTJ that f° r large times, r ~ 1,000, 
the three models converge to the same power law decay as was anticipated by 
Eqs. (I43p . (I5TT) . and (1561) with a + 1 = 1.94, which differs from the behavior 
of the empirical pausing-time density shown in Fig. [TTJ This is not surprising 
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Fig. 11. The pausing time distributions for the U.S. dollar /Deutsche mark. We 
compare the pausing time distributions obtained from the Laplace inversion of 
Eqs. (|4ip . (|49p and (|54p with the empirical pausing time distribution (|34p . The 
results produced by the three ansatzes are polotted with the value of a obtained 
from the fit of fewer days in Fig. [10] with r given by Eq. (|42p . The empirical fit takes 
r = t/{T). Note that both normalized times are quite similar. The inset shows in a 
logscale how the three models collapse in the same power law when r > 1,000. 

since this time scale is outside the region of validity of the empirical distri- 
bution (see Fig. [5] and Eq. Indeed, if we assume that t ~ t(T), then 
r = 1, 000 approximately corresponds to seven hours of trading, an interval of 
time exceeding the trading day. 



5 Conclusions 

We have described some uses of CTRW formalism in the context of finan- 
cial markets. The CTRW formalism provides insight that is able to relate not 
only the market microstructure activity with the probability distribution of 
intraday prices, but also with the distributions of daily, weekly, or longer- 
time prices. Techniques based on the CTRW also can be applied to inverse 
problems in which it is possible to analyze the microstructure of financial pro- 
cesses, represented by their transaction statistics, based only on low-frequency 
(daily, weekly, etc.) data. In this way, we can have some knowledge of the mar- 
ket microstructure by taking the information related to low frequency price 
changes. 
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The formalism, in its direct version (i.e., from high-frequency data to lower 
frequency statistics), depends on estimates of the pausing-time density, ip(t), 
and the jump density, h(x). A first and simpler approach consists in assuming 
that these densities are independent. However, in many practical situations 
one certainly expects some degree of correlation, in the sense that large return 
increments are infrequent. We have proposed two possible ways to model this 
correlation. We have described some general features of the formalism that 
hold for any model of correlation between sojourn times and return jumps. 
Thus (i) the variance grows linearly with time as t increases, and (ii) the tails 
of the return distribution p(x, t) have the same appearance than those of the 
jump density h(x). We have also applied the direct formalism to data provided 
by the U.S. dollar/Deutsche mark future exchange rate. High frequency data 
indicate that ip(t) and h(x) are well described by power- law densities. We have 
therefore shown that the volatility has a diffusion-like behavior at long times 
and that the tails of the return distribution, p(x,t), follow a power law with 
the same exponent as that of h(x). 

We have applied the inverse formalism to daily data provided by four compa- 
nies: IBM, Gillette, Dupont and 3M. We have focused on the volatility and 
its behavior as a function of the time lag. The data for all companies exhibit 
a slight but clearly anomalous diffusive behavior at short times since the vari- 
ance can be fitted with a polynomial function growing as t a in the short time 
regime (few days), and as was expected, the variance describes an ordinary 
diffusion-like behavior for return changes of a week or more, that is, a linear 
dependence on time for large time lags. All of this implies that the shortest- 
day behavior of the unknown pausing-time density ip(t) follows a power law 
consistent with high frequency data for the US-Deutsche mark case and that 
for larger times the pausing-time density decays exponentially. 

Finally, these results seem to indicate the suitability of the CTRW framework 
to study other complex issues in finance. We can mention that the present 
formalism can be useful in the study of problems related to large fluctuations 
and extreme values central to risk control and also to exotic derivative pricing. 
All of these are under present investigation. 
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A Asymptotic results 



Let us prove the results given in Section 13.11 Since AX and T are random 
variables representing respectively returns increments and sojourn times, we 
can write p(u, s) in the form (see Eq. (J3])) 

p(u, s) = (exp{iuAX — sT}). 

Expanding this average around s = 0, we get 

p(cu, s) ~ h(u) - s(Te iujAX ) (s < (T)- 1 ) , (A.l) 

where 

h(u) = ( e iujAX ). 

On the other hand, if we assume that the sojourn-time density ip(t) has a 
finite first moment, then 

$(s) ~ 1 - s(T) (s « (T)- 1 ) , (A.2) 

where (T) is the mean sojourn time. Therefore, from Eq. ffTUl) we have 

p (u,s)^(T) ( S «(T)- 1 ). (A.3) 



Substituting Eqs. flATB and (TAT3]) into Eq. © yields 



P(uj,s)~ = ^ (s^iT)- 1 ). (A.4) 



By virtue of Tauberian theorems (Handelsman and Lew 1974), the asymptotic 
behavior of the characteristic function p(u, t) as t — » oo will be given by the 
inverse Laplace transform of Eq. (1A.4j) . The resulting inverse reads 



We emphasize that the asymptotic expression in Eq. flA.5[) is totally general 
regardless of the correlation between sojourn times T and return increments 
AX, requiring only that the mean sojourn time (T) be finite. Obviously for a 
usable form of this asymptotic expression we will need to know the correlation 
between T and AX. In other words, we will have to specify a form for p(x, t) 
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since the correlation appearing in Eq. (1A.5I) is related to p by 



(Te^ AX ) = - |-p(o;, s) . 

s=0 

Thus, for instance, making use of ansatz in Eq. ( fT4"l) . we have 

(Te wAX ) = (T)h 2 (u) 



so that 



p(uj,t) 



1 



exp < ~ — 

h?{u) l \ (T)h?{uj) . 



(t » (T» 



(A.6) 



Now assume that /i(a;) has a finite second moment H2 so that /i(cj) ~ 1— /i 2 ^ 2 /2 
as yU2^ 2 "C 1. In this case, for values of uj sufficiently small, Eq. (IA.5I) yields 
the Gaussian density 

which is equivalent to the result given by the central limit theorem. 

On the other hand, for a long-tailed jump density for which h{u) ~ 1— k\oo\ a for 
small to (note that if < a < 2, h(x) has infinite variance) and assuming that 
for uj small (Te tujAX ) ~ (T), we obtain from Eq. (IA.5I) the Levy distribution 



We finally show that the tails of the return pdf p(x, t) are the same as that of 
the jump pdf h(x). Indeed, as is well known the tails of p(x,t) as \x\ — > 00 
are determined by the behavior of its characteristic function as uo — ► (Weiss 
1994). Thus suppose that uj is small so that (Te iujAX ) ~ (T) and t is moderate 
in the sense that the approximation given in Eq. flA.51) is still valid but in such 
a way that [1 — h{ui)]t/ (T) 1. Then from Eq. flA.51) we see that for small u, 



p{(j,t)~-[l-h(u)]t/(T), 
which for large values of \x\ is equivalent to 

p(x,t) ~ -±-h(x), (t«<T». 

Therefore, we have proved Eq. (124")) in item (c) of Section [3TT1 and shown that 
the tails of p(x,t) are equal to the tails of h(x). 
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B Pausing time density limits 



In order to get the asymptotic behavior of ip(t) as t — > oo and as t — > 0, we first 
write the Kummer function appearing in Eq. (I4"9"j) in terms of the incomplete 
gamma function: U(a,a,x) = e x T(l — a,x) (Magnus et al. 1966). Therefore, 

m = [kT ^J /a e ^)/^ v {l _ 1/a> fcr ( a)/Al3S «] . (B.l) 
M s 



Now, using Tauberian theorems (Handelsman and Lew 1974), the behavior 
as t — > of ijj(t) will be given by the behavior of ip(s) as s — > oo while the 
behavior of ip(t) as t — > oo is determined by V'C 5 ) as s — > 0. 

(a) Short-time behavior. We use the following series representation of the in- 
complete gamma function: 
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where 

/3 = 1 - l/a, x = kT(a)/fis a , 
and Lg) (x) are the generalized Laguerre polynomials (Magnus et al.). Since 
L(f)(0) = T(l +/3 + n)/n!r(l + /?), we have " 
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Ac" + O (x^ 1 ) 



where 
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(B.2) 



Hence 
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which implies (Handelsman and Lew) 
Ak 



A*2 



(B.3) 



(b) Long-time behavior. In order to obtain the long-time behavior of the den- 
sity ip(t) corresponding to the volatility, we employ the following asymptotic 
behavior for the gamma function (Magnus et al.): 



T(P,x) =e~ x x p ~ 
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Thus 

and (Handelsman and Lew) 

m ~ — f 2 r 1 -" (r » ^) . 

ajr(l + a) 
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